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Abstract. The production of dijets with a leading neutron in ep-interactions at HERA is calculated in 
leading order and next-to-leading order of perturbative QCD using a pion-exchange model. Differential 
cross sections for deep-inelastic scattering (DIS) and photoproduction are presented as a function of several 
kinematic variables. By comparing the theoretical predictions for DIS dijets to recent HI data, the pion flux 
factor together with the parton distribution functions of the pion is determined. The dijet cross sections in 
photoproduction show factorization breaking if compared to the HI photoproduction data. The suppression 
factor is S = 0.48 (0.64) for resolved (global) suppression. 



PACS. 12.38.Bx Perturbative QCD calculations - 13.60.-r Photon interactions with hadrons 

1 Introduction 

In recent years, the validity of QCD factorization in hard diffractive scattering has become an important issue experi- 
mentally and theoretically. Factorization in hard diffraction means that the observed cross sections in hard diffractive 
processes can be calculated by a convolution of diffractive parton distributions with parton- level cross sections. Hard- 
scattering factorization has been proven by Collins 1 for inclusive diffractive deep- inelastic scattering (DIS), i.e. for 
the diffractive structure functions. It is supposed to be valid also for subprocesses like jet production and heavy-quark 
production in the DIS region. The proof of the factorization formula also appears to be valid for the direct photo- 
production of jets and heavy quarks 1 . Factorization does not hold for hard processes in diffractive hadron-hadron 
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scattering. Here, soft interactions between the two hadrons, and their remnants, occur in both the initial and final 
state, which prevents using the same steps as in the proof for diffractive DIS. Therefore, factorization fails also for 
resolved photoproduction. The failure of factorization in hadron-hadron scattering is observed experimentally [2] . The 
cross section for diffractive dijet production at CDF is suppressed relative to the prediction based on diffractive parton 
distribution functions (PDFs) from the HI collaboration [3] by one order of magnitude |2J- The breaking of hard factor- 
ization in diffractive dijet photoproduction is also firmly established by analyses of HI 4 and ZEUS [5] experimental 
data. Whether these experimental data are consistent with the breaking in the resolved component alone or whether 
this breaking occurs also in the direct photoproduction cross section is still not satisfactorily proven. It seems that the 
data are better described by a global suppression of the direct and resolved contribution by about a factor of two. It 
is important to note that this suppression is only visible if the data are compared to the next-to- leading order (NLO) 
QCD predictions as it was first shown by us in [5][7]. 

Factorization breaking is expected not only in the diffractive region, xp <C 1, where xp is the momentum fraction 
transferred to the exchanged particle in the t-channel, but also at larger values of xp, where Regge exchanges other 
than the pomeron occur between the initial and final state proton. These secondary Regge-pole exchanges are not 
only present for p —> p transitions but also in p — > n transitions, in which pion exchange is strong. Therefore dijet 
photoproduction with a leading neutron could also be a candidate for factorization breaking as was already suggested 
in 6 . Dijet photoproduction, e + + p — > e + + n + jet + jet + X' ', with pion exchange has been studied in leading order 
(LO) and NLO in [5j and compared to ZEUS experimental data UJ. Recently the HI collaboration has measured 
these dijet cross sections for photoproduction (Q 2 < 0.01 GeV 2 ) and DIS (Q 2 > 2 GeV 2 ) The photoproduction 
cross sections were compared to the NLO predictions of |SJ, and good agreement was found concerning the shape and 
normalization of distributions for various kinematic variables similar to the comparison done in [SJ to the ZEUS data. 
Neither in [S] nor in |l(Jj the factorization breaking of the dijet photoproduction cross sections has been investigated. 
This is difficult, since the breaking shows up dominantly in the normalization of the cross sections and to a much 
lesser extent in the shape of the distributions. 

In the pion-exchange model, the normalization of the neutron-tagged cross sections depends first on the splitting 
function of a proton into a pion and a neutron fTr/p(xL,t). Here, xl and t are the two variables which describe the 
proton- neutron vertex, xl is the fraction of the initial-state proton energy transferred to the neutron, and t is the 
square of the momentum transfer between the proton and neutron. Second, the normalization depends on the parton 
distribution functions of the pion, for which several models exist in the literature 11 . 
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The pion flux fn/p(zL,t) can, in principle, be measured in charge-exchange processes in soft hadronic reactions, 
where an initial-state proton is transformed into a final-state neutron, p — ► n, with small momentum transfer. A 
successful phenomenological description of the corresponding data has been given in the framework of reggeized 
isovector exchanges, such as 7r, p and a-i with the pion dominating the p — > n transition, in particular at small squared 
momentum transfer t between proton and neutron |12j . Unfortunately these cross sections are not given by pure Regge- 
exchange amplitudes. They are modified by soft rescattering of the incoming and outgoing hadrons, which influences 
the normalization of the cross sections, i.e. it leads to modified flux factors f n /p(xL, t) due to the absorption of leading 
neutrons and the ingoing protons, which depends on xl and t. Such corrections were studied some time ago by several 
authors |13| . 

The pion PDFs are constrained by di-muon and prompt-photon production data from fixed target experiments 
that are sensitive to the valence quark distribution in a Bjorken-x range relevant for dijet production on pions 
Unfortunately, these constraints are not very restrictive so that the existing parameterizations jl lj differ in the relevant 
x-range, which leads to appreciable differences in the calculated dijet photoproduction cross sections 

The most direct way to determine the flux factor times the structure function of the pion F$(x, Q 2 ) is to measure 
the inclusive DIS cross section with a tagged neutron. On the basis of Collins' factorization proof, we expect in this 
case no absorptive interactions. Such cross sections for leading neutrons have been measured in ep-collisions at HERA 
[T%| . Unfortunately these data have not been analyzed towards determining the PDFs of the pion, assuming a fixed 
ansatz for the p — > n vertex, as it has been done towards the PDFs of the pomeron from diffractive semi-inclusive cross 
sections |51ll5| . Therefore we shall follow a different route. Assuming one of the PDFs of the pion in the literature, 
for example those of Gliick, Reya and Vogt (GRV) we determine the normalization of f v / p {xL,t) by comparing 
the NLO dijet cross sections for Q 2 > 2 GeV 2 to the data of the HI collaboration ^U|. With this f n /p{xL,t) in 
conjunction with GRV's pion PDFs, we have calculated the NLO dijet cross section for photoproduction (Q 2 ~ 0) on 
the basis of |S] and compare them with the experimental results in jl(J) . Here the experimental tagged neutron DIS 
and photoproduction cross sections are obtained in a common analysis in the same region of the neutron kinematic 
variables. If we assume that the dijet production data in the DIS region are not influenced (or influenced very little) 
by absorptive corrections, i.e. have no (or only a very small) breaking of factorization, then the failure of our NLO 
photoproduction cross sections to describe the corresponding data of 10. will give us the amount of factorization 
breaking for photoproduction. 

In Section 2, we shall describe the calculation of the dijet cross sections with a leading neutron, together with the 
kinematic variables, and define our input for the pion flux and the pion PDFs. The NLO cross sections are compared 
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Fig. 1. Generic Feynman diagrams for the scattering process e + p — *• e + n + X (left) and for the production of two jets in the 
one-pion exchange model (right). 

to several measured DIS cross sections from JH]; so that our assumptions concerning f n / p (xL,t) and the pion PDFs 
can be tested. Section 3 contains our results for the dijet cross sections in the Q 2 ~ region and the comparison with 
the experimental data of ^UJ- O n this basis we test also whether the factorization breaking can be attributed to the 
resolved component alone. In Section 4, we give a short summary and draw some conclusions. 



2 Dijet Cross Sections in DIS 
2.1 Kinematic Variables and Input 

The event kinematics has already been described in QJjj. Here we recall the definition of those variables which are 
needed in the calculations of the cross sections. The reaction 

e+(k) + p(P) -» e+(k') + n(P') + X(p x ), (1) 

where X is the hadronic system containing at least two jets, is characterized by the four-momenta k and k' of the initial 
and scattered positron and by P and P' , the four- momenta of the ingoing proton and outgoing neutron, respectively, 
as sketched in Fig. 1 (left). The positron-photon vertex is described by the exchanged photon virtuality Q 2 and the 
positron's inelasticity y, 

Q 2 = -Q 2 = -(* - k'f and y = ||. (2) 
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In the HI experiment 10 1, 2 < Q 2 < 80 GeV 2 and 0.1 < y < 0.7 in the DIS region and Q 2 < 10~ 2 GeV 2 and 
0.3 < y < 0.65 for the photoproduction selection. The protons at HERA have the energy E p — 820 GeV and collide 
with E e = 27.6 GeV positrons. 

The two variables, which describe the proton- neutron vertex, are the fraction Xl of the energy of the initial-state 
proton E p carried by the neutron and the square of the momentum transfer t between the proton and the produced 
neutron, 

P'q E n 

if, = ~ and 

Pq Ep 

_ ~ P T n ,, , / »',, 2 



t=(P-P'y~-^-{l-x L )[-i-m}), (3) 

XL V X L J 

where E n is the neutron energy, px n is the momentum component of the neutron transverse to the direction of the 
ingoing proton, and m n and m p are the neutron and proton masses, respectively, xl and t are determined from the 
measured energy and scattering angle of the leading neutron. In the HI experiment E n > 500 GeV and 6 n < 0.8 mrad. 

In the pion-exchange model, the photon interacts with a pion emitted from the proton. In this model, the process 
e + +p — > e + + n + jet + jet + X' as sketched in Fig. 1 (right) is described by the variable which, neglecting masses, 
is the fraction of the four- momentum of the pion q' = P — P' participating in the hard interaction. It is related to x p , 
the fraction of the four-momentum of the proton, which enters the hard interaction, by x p = — Xl). ij, is related 
to Xp introduced in Section 1 in connection with diffractive jet production via xl = 1 — xp. 

In the scattering of 2 — > 2 massless partons, the fractions of the four-momenta q = k — k' and q' = P — P' transferred 
to the partons are given by 

Jet = (4) 



and 

™Jet 

2^(1 - x L ) ' 



so that 

J- = U^re _ (g) 



The sums in Eqs. (4)-(6) run over the variables of the two jets in the final state. and if et denote the transverse 
energies and pseudorapidities in the laboratory system. Strictly speaking, Eq. (4) is correct only for photoproduction, 
where q 2 ~ 0, and furthermore q 12 = 0. The energy fraction contributing by the exchanged virtual photon to the 
production of the dijets is x J 7 ct , whereas the corresponding contribution of the virtual pion (or possibly of a reggeized 
p or 02) is x^ 1 . In Eq. (4), E 1 = yE e is the energy of the ingoing virtual photon. In NLO, also three jets can be 
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produced in the final state, but these contributions have been removed from the theoretical prediction in accordance 
with the experimental analysis, which contains only an exclusive dijet sample (see below). 

As in the jet analysis of the experimental data we use the cone algorithm with radius R' = 1 JHI f° r 
the jet definition and the combination of two partons into one jet in the NLO contributions of the DIS and the 
photoproduction sample. In our previous work [8], we had used the fcr-cluster algorithm instead. The jet finding is 
performed in the j*p center-of-mass frame with transverse energies defined relative to the 7* momentum axis. The 
laboratory pseudorapidity of each jet is restricted to the region — 1 < ?/ Jct < 2. The transverse energies of the two 
jets with the largest Et are constrained to the region > 7 GeV and E^ t2 > 6 GcV for both the DIS and the 

photoproduced jets in accordance with the requirement of asymmetric Ej, '-cuts |17| . The calculated cross sections are 
restricted to the selection of exactly two jets, i.e. the rather small contribution of three jets with E^ t3 > 6 GeV is not 
included as in the experimental selection jl(Jj . 

For the production of dijets in the DIS region, we consider only the contribution of the directly coupled 7*, although 
at Q 2 as low as 2 GeV 2 the resolved contribution might be relevant U^]. For the photoproduction of dijets we take 
into account the direct and the resolved process both in LO and NLO. The observable x J 7 ot is sensitive to the amount 
of direct and resolved processes. 

In the pion-exchange model, the cross section for j*p scattering to the final state nX (see Eq. (1)) takes the form 

d<j(l*P -> nX) = U /p (x L ,t) da( 7 *7T+ -> V), (7) 

where f^/p is the pion flux for the transition p — > n + tt + and d(r(7*7r + — > V) stands for the cross section of the hard 
7*-7r + interaction. The splitting function or pion flux is usually parameterized by different forms. We choose 

M*r.,t) = ^7-^(1 -*Lr^[F{ XL ,t)f. (8) 

Here, g n -wp is the coupling constant of the nirp vertex, is the pion mass and a„-(t) = a'(t — m 2 .) is the Regge 
trajectory of the pion. F(xL,t) is a form factor, which describes the off-shell behavior of the virtual pion and/or 
possible final-state rescattering effects of the neutron. We choose the so-called light-cone form factor 

F{x L ,t) = exp[R 2 (t - ml) /{I - x L )}. (9) 

This choice is usually associated with the flux without Regge trajectory factor, i.e. a' = in Eq. (8) (T^|. The pion- 
nucleon coupling constant is known from low-energy 7rV and NN scattering data. We take g^p/Air = 2 x 14.11 |2(J| . 
Other choices of F(xl, t) have been used in the literature. We consider only the form in Eq. (9), since it has been used 
also in [TU|. 
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Another important input are the PDFs of the pion. For this several choices are available in the literature We 
choose the parameterization of Gliick, Reya and Vogt (GRV) which was also chosen in ^Uj. They provide PDFs 
in LO and NLO and also include the charm contribution. For the calculation of the photoproduction cross sections we 
need the PDFs of the photon for the resolved part. A popular parameterization is the one of Gliick, Reya and Vogt 
(GRV) which was also the choice in ^Uj. The Aj^g parameter, which we need in the LO and NLO formulas for 
a s , is adjusted to the PDG 2006 edition |22] value for a s {m z ) = 0.1176. This yields for four flavors A^L = 0.118 
GeV in LO and A^ = 0.307 GeV in NLO. Of course, our NLO predictions for DIS and photoproduction dijet cross 
sections depend on this value. For the LO (NLO) predictions we use LO (NLO) hard scattering matrix elements with 
the one-loop (two-loop) formula for a s and the Aj^g values given above. Unfortunately, the /Lvalues used for the 
evolution of the GRV pion and photon PDFs are somewhat different, namely 0.200 GeV in LO and NLO. In the 
case of photoproduction, the photon flux is calculated with the usual Weizsacker- Williams approximation including 
the non-logarithmic corrections as calculated in 23;. The renormalization and factorization scales are equal to the 
maximum transverse energy of the outgoing jets. Now, all parameters and PDFs are specified except for the radius 
R in the light-cone form factor in Eq. (9) . The value of this parameter will be fixed by comparison of the theoretical 
cross section in DIS dijet production with the data of ^Ujj i- e - with the measured cross sections da/dE^ , da/drf et , 
da / dlog 10 (x^ t ) , and da/dQ 2 . In |T|j], the photoproduction dijet data have been described very well with the choice 
R = 0.65 GeV^ 1 . 

In addition to pion exchange, secondary Regge exchanges (as for example p and 02) are possible. Such contributions 
could be disentangled by a careful study of the flux factor as a function of xl and t (or px n ) ■ Such data are not available. 
In |10|. the neutron energy dependence of the DIS and photoproduction dijet events were measured. The uncorrected 
data were compared to a pion-exchange model within a Monte Carlo simulation that included detector effects. The 
photoproduction data are reasonably described in shape and magnitude, whereas the DIS data are reproduced quite well 
in shape, but somewhat overestimated in magnitude. We take these results for a good indication that pion-exchange is 
dominating over p- and (^-exchanges at least in the region of xl > 0.6. Since additional p- and a2-exchanges influence 
only the p — > n flux factor, this factor should be the same for the DIS and photoproduced dijets due to factorization. 

2.2 DIS Dijet Cross Sections 

The calculations of the dijet cross sections in the DIS region have been performed with the NLO Monte Carlo program 
JETVIP 24 . This program calculates jet cross sections in LO and NLO in DIS using the so-called phase space 
slicing method with an invariant mass cut-off to cancel the infrared and collinear singularities. The program, originally 
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constructed for jet production in 7*p interactions, could easily be modified for our purpose to calculate jet cross 
sections with a tagged neutron. The only change is that the PDFs of the proton are replaced by the PDFs of the pion 
times the pion flux. 

Since the neutron kinematics is not fixed in detail by the experiment, we had to integrate over a finite region in 
xl and t according to Eqs. (3). We did this in accordance with the specifications of the HI experimental analysis. 
Except for the outgoing positron and the leading neutron, the final state consists of two or three jets in addition to the 
remnant jet of the pion. The two-jet sample contains the bare parton jets from the LO and virtual NLO contributions 
and the two jets originating from the recombination of two partons in the three-parton contribution using the cone 
algorithm. Then we calculated the differential cross section d 2 a/dET<ir]. where Et and rj are the transverse energy 
and the rapidity of the jets in the two- or three-jet sample with E^ 1 > 7 GeV and E^ t2 > 6 GeV and E^ t3 < 6 GeV, 
i.e. from the three-jet sample the events with a hard jet with E^ t3 > 6 GeV are left out. This defines the exclusive 
dijet sample as in the analysis by HI ^Uj. The cross section da j AEt^j] was integrated over the region — f < rf ct1 ' 2 < 2 
to yield da/dEx and integrated over Et with the asymmetric Ex-cut defined above to give us da/dif et and similarly 
for dcr/dlog 10 (4 et ) and da/dQ 2 . 

The results for da/dE 1 ^" , der/d^ ', dcr/dlog 10 (x J 7 f !t ) and da/dQ 2 are shown in Figs. 2a-d. In these figures we have 
plotted the experimental data from and three theoretical predictions, in LO and NLO with R = 0.65 GeV -1 
and one in NLO with R — 0.55 GeV -1 . For R = 0.55 GeV -1 we show also the scale variation in NLO by varying 
the scale in the interval fi/2 to 2fi. It is seen quite clearly that the NLO prediction with R = 0.65 GeV -1 does not 
describe the data, whereas the prediction with R = 0.55 GeV -1 is in reasonable agreement with the measured cross 
sections. We can quantify these statements by calculating the Xdf f° r the two cases. For R = 0.65 (0.55) GeV -1 
we get Xdf = 1-78 (0.78) for 24 degrees of freedom. We observe that the NLO corrections are very important in all 
four distributions, as one can see by comparing the predictions for R = 0.65 GeV -1 in LO and NLO. Our theoretical 
predictions will be modified by hadronic corrections, which stand for the difference between jets from hadrons, on which 
the measured cross sections rely, and jets built out of partons which are the subject of the theoretical predictions. 
Unfortunately these hadronic corrections, which are usually obtained from Monte Carlo programs, which simulate LO 
cross sections with parton showering and parton-to-hadron transitions, are not available for the DIS dijet cross sections 
with tagged neutrons. From experience with such cross sections with tagged protons, we know that these corrections 
are of the order of 10 — 20% except for the x J 7 et -distribution, which is very much changed by these corrections. Actually, 
da/dxj has also been measured in |10) for DIS dijets. Due to the sensitivity of this distribution to hadronic corrections, 
we did not consider this cross section suitable to give us the a trustworthy i?-value. 
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Fig. 2. Differential cross sections for deep-inelastic dijet production with a leading neutron. The HI data are compared to 
perturbative QCD predictions in LO and NLO and with two different values for the pion-neutron radius R. 



3 Photoproduction Cross Sections 



As in , the calculation of the photoproduction cross sections is based on the formalism fully described in our previous 
work [53]. The cross sections which we shall evaluate are the same as in the DIS case described above. However, we 
now include also dcr/da; J 7 ct in order to see the contribution of the resolved part in the region x J 7 ct < 1 more clearly. In 
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Fig. 3. Differential cross sections for dijet photoproduction with a leading neutron. The HI data are compared to perturbative 
QCD predictions in LO and NLO and with two different values for the pion-neutron radius R. 



the following we show the results for da/dE^ 1 , do-/dif ct , da/dx^ and d<7/dlog 10 (a;|f :t ) in Figs. 3a-d. In these figures, 
four different predictions are plotted and compared with the data from |10|. Three predictions are in NLO, one with 
the radius R = 0.55 GeV -1 , where hadronic corrections |ll)U2t>) are included, and two with R — 0.65 GeV -1 with 
and without hadronic corrections, and one is in LO with R = 0.65 GeV -1 without hadronic corrections. The LO 
prediction is far off the experimental data showing that only NLO predictions are relevant and that the X-factor (K = 
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NLO/LO) is large. For the Eji '-distribution in Fig. 3a, the prediction with R = 0.55 GeV -1 , where the shaded band 
gives the scale-dependence of the cross sections, lies always above the data points except for the two largest Ex-hms. 
The prediction with R = 0.65 GeV -1 agrees much better with the data in agreement with the results of ^UJ- I n (QS 
the cross sections are somewhat smaller, since there a s is smaller due to the choice of a smaller A^— value (0.200 

' s MS y 

GeV). From this figure, it is clear already that the radius R = 0.55 GeV -1 gives a bad description of the HI data and 
a reasonable agreement would be possible only for R > 0.65 GeV -1 . This is seen even more clearly in Figs. 3b,c,d, 
where the if et , and log 10 (x^ t ) distributions are compared to the measured cross sections in JU]. The predictions 
lie above the data points, even if the theoretical error estimated by the scale dependence is taken into account. To 
reproduce the data in Figs. 3a-d a radius larger than R — 0.65 GeV" 1 would be needed. However, the radius R = 0.55 
GeV -1 was fixed by the DIS dijet cross sections. Our results in Figs. 3a-d, when compared to the photoproduction 
data of the HI collaboration, demonstrate that they can not be explained with the same pion flux plus pion PDFs. 
i.e. with R — 0.55 GeV" 1 . This shows us that factorization breaking is present in photoproduction of dijets with a 
leading neutron with a breaking factor of S ~ 0.6. This breaking factor changes somewhat from cross section to cross 
section and as a function of the kinematic variables. 

The details are shown in Figs. 4a-d, where we show the ratio HI Data/Theory for the four cross sections with 
R = 0.55 GeV -1 , with and without hadronic corrections. In addition to the ratios following from the results in Figs. 
3a-d, we show in Figs. 4a-d also another ratio, where the theoretical prediction is again obtained with R = 0.55 
GeV -1 including hadronic corrections, but now with a suppression factor S = 0.48 applied to the resolved cross 
section together with the /x 7 scale-dependent part of the NLO corrections to the direct cross section. Here, /x 7 is 
the factorization scale of the photon vertex. This additional suppression in the direct cross section introduced in |27j 
serves the purpose to eliminate the /z 7 -dependence of the sum of the direct and resolved cross section. We see from 
Figs. 4a-d that this ratio lies near to one if we take into account the scale variation of the theoretical prediction. The 
^'-distribution yields values for this ratio larger than one (except for the last bin), which overlap, however, with the 
scale variation band. It is known that the ^'-distribution suffers from rather large hadronic corrections as seen for 
example Fig. 3c, where the correction in the bin a; J 7 et € [0.6; 0.8] is as large as a factor of two. 

The two suppression factors, S ~ 0.6 for the full direct and resolved contributions or S = 0.48 for the resolved 
and initial-state singular part of the direct contribution, can be compared with the suppression factor obtained in 
|28| . In this work, the spectra of leading neutrons, both in photoproduction and in DIS, were studied and compared 
to recent ZEUS experimental data (Tl|. It was found that the photoproduction cross section on the basis of the 
pion-exchange model agreed with the data if it was reduced by a factor of about 0.4. This factor could be quite well 
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Fig. 4. Ratios of HI Data over NLO QCD for dijet photoproduction with a leading neutron, with and without hadronization 
corrections and including (full points/lines) a suppression factor S for the resolved component and its scale-dependent direct 
NLO counterpart. 
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explained by absorptive corrections to pion exchange and can be compared with the suppression factor of 0.48 for 
the resolved/direct-IS components. Furthermore, it must be emphasized that the absorptive corrections for the total 
photoproduction cross section with a tagged neutron may differ from the absorptive corrections for the production of 
a pair of high-Er jets, even in the resolved case. 

In Fig. 5, we show separately the Ej. '-distribution as measured by the HI collaboration and compare it to theoretical 
NLO QCD predictions including hadronization corrections and using the radius R = 0.55 GeV -1 . In the full histo gram, 
we test the hypothesis of suppressing only the resolved (and direct initial-state) contributions with a factor S = 0.48 
and obtain quite a reasonable value of Xdf = 1-2- The shaded band indicates again the theoretical uncertainty coming 
from simultaneous variation of the renormalization and factorization scales by a factor of two around the central scale. 
It overlaps with the experimental error bars in all but the highest two bins, whose precision is obviously limited by 
statistics. In contrast, the dashed histogram tests the hypothesis of suppressing all direct and resolved contributions 
with a fitted factor S = 0.64, as proposed in the HI analysis of diffractive dijet photoproduction, leading to a higher 
value of Xdf ~ l-7< We take this as an indication that the first hypothesis describes this most robust (and exponentially 
falling) distribution better. 

Finally, we wish to comment on our earlier analysis |Hj of the photoproduction dijet cross sections as measured 
by the ZEUS collaboration several years ago At this time, data on dijets with tagged neutrons in DIS were not 
available. Therefore, we fitted the pion-neutron radius R of the light-cone form factor to the ZEUS data with the 
result R = 0.5 GeV -1 . As shown in this work, the HI photoproduction data JU| are badly described with this radius. 
The reason for this mismatch is the fact that in [5] we have chosen different parameterizations of the pion PDFs, 
i.e. SMRS versus GRV 11 , and of the photon PDFs, i.e. GS96 [H| versus GRV [2T]. In particular, the choice of the 
pion PDFs has a strong influence on the absolute value of the cross section (see for example Fig. 3 in |S]). The ZEUS 
collaboration demonstrated in 2002 JT1J that the shape of their measured F% distribution agrees quite well with the 
GRV parameterization and much less with SMRS. Furthermore, the normalization of the theoretical dijet cross section 
was influenced by the choice of the difference in the £r-cuts for the two jets, which was set to zero in the analysis of 
the experimental data. 

4 Conclusion 

In Summary, we have performed a comprehensive NLO QCD analyis of dijet production with a leading neutron in both 
deep-inelastic scattering and photoproduction. We emphasized the question whether factorization breaking occurs not 
only in diffractive photoproduction, but also in photoproduction with p — > n transitions. Assuming that the latter 
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Fig. 5. Ej, -distribution of photoproduced dijets with a leading neutron. The HI data are compared to NLO QCD calculations 
corrected for hadronization effects in the pion-exchange model with pion-neutron radius R = 0.55 GeV -1 and suppression of 
the resolved/direct initial-state contributions (full) and of all direct and resolved photon contributions (dashed). 
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are dominated by one-pion exchange and well described by GRV pion PDFs and a light-cone flux factor, we were 
able to fit its free parameter, the pion-neutron radius R, to recent HI DIS data with the result R = 0.55 GeV" 1 for 
Xdf = 0-78 an d 24 degrees of freedom. In contrast, a radius of R = 0.65 GeV -1 as used by the HI collaboration leads 
to a considerably larger Xdf °f 1-78. 

When applying the fitted radius of R = 0.55 GeV -1 to our NLO QCD predictions for dijet photoproduction with 
a leading neutron, we seriously overestimate the corresponding HI data. We therefore conclude that factorization 
breaking occurs not only in diffractive ep — > e'pX, but also in ep — > e'nX scattering. Only after including absorptive 
corrections, a good phenomenological description of the HI data can be obtained. A suppression of resolved photopro- 
duction and its factorization-scale dependent direct NLO counterpart only with a suppression factor of S — 0.48 seems 
to be favored over a global suppression of all direct and resolved contributions with a suppression factor of S = 0.64, 
at least in the steeply falling i/^'-distribution, where we obtain values of Xdf = 1-2 versus 1.7. 

M.K. gratefully acknowledges the hospitality of the Institute of Theoretical Physics at the University of Gottingen, where part 
of this work has been completed. 
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